Dynamics of a Heisenberg spin chain in the quantum critical regime: NMR 
experiment versus effective field theory 
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A comprehensive comparison between the magnetic field- and temperature-dependent low fre- 
quency spin dynamics in the antiferromagnetic spin-1/2 Heisenberg chain (AFHC) system copper 
pyrazine dinitrate, probed via the 13 C-nuclear magnetic resonance (NMR) relaxation rate T-j -1 , and 
the field theoretical approach in the Luttinger liquid (LL) regime has been performed. We have 
found a very good agreement between the experiment and theory in the investigated temperature 
and field range. Our results demonstrate how strongly the quantum critical point affects the spin 
dynamics of Heisenberg spin chain compounds. 

PACS numbers: 75.10.Pq, 71.10.Pm, 76.60.-k 



The occurrence of quantum phase transitions (QPT) 
in systems of correlated electrons is a very important 
topic in current solid-state physics. These transitions are 
present in, e.g., high-T c superconductors, heavy- fermion 
metals, or magnetic insulators [l|. The phase diagrams 
of systems from the first two classes are in general com- 
plex due to several interaction mechanisms. In contrast, 
the purely magnetic interactions in magnetic insulators, 
in particular for one-dimensional spin systems, give the 
rare occasion to perform exact calculations of their char- 
acteristics and to compare them with experimental data 
sets of well characterized sample systems |2|]. The spin 
properties of organic-based low-dimensional magnets can 
be fine-tuned by chemical synthesization methods. This 
well controlled synthesis allows the systematic investi- 
gation of the magnetic properties with well established 
methods such as neutron scattering, ESR, dc/ac magne- 
tometry, /iSR or NMR. In low-dimensional magnets the 
high sensitivity of NMR |3| to local hyperfine fields allows 
to perform detailed studies of, e.g., phase transitions, or 
the local distribution of spin moments. 

The isotropic AFHC model is one of the main 
paradigms of quantum many-body physics both from the 
experimental and theoretical viewpoint. Its static charac- 
teristics were successfully compared with experimentally 
studied features of quasi-one-dimensional magnetic com- 
pounds, synthesized recently Q. For dynamical prop- 
erties, especially in the vicinity of the quantum critical 
point, there is still an insufficiency of experimental data 
sets for the low frequency part and their comparison to 
calculations. 

In this Letter we present a detailed comparison of field 
theory results with the data of recent NMR experiments 
|5| probing the spin dynamics in a wide field- and temper- 
ature range in one of the best realizations of the AFHC 



model, namely Cu(C4H4N2)(N03)2 (also known as cop- 
per pyrazine dinitrate or CuPzN) [fj . We find an extraor- 
dinary good agreement in the low temperature behavior 
across the field driven QPT. 

Due to a relatively low value of the coupling constant 
J/k B =W.7 K in CuPzN, the critical field value B s = 14.6 
T is well reachable by standard laboratory equipment. 
Therefore, one can examine the spin dynamics in the re- 
gion of fields and temperatures, where spin-spin corre- 
lations manifest themselves in the most prominent way, 
and compare with the results of various theoretical meth- 
ods. The inter-chain interactions are supposed to be 
small, so that the magnetic ordering (T c ~ 107 mK) did 
not affect the AFHC behavior down to the lowest T stud- 
ied in the NMR-cxpcrimcnt. Whereas these data were 
compared with numerical Quantum Monte Carlo (QMC) 
simulations, their agreement with results of a field the- 
ory approach (which serves as a very good description 
namely at low T, where QMC simulations often produce 
larger errors) for the LL regime is checked in our work. 
A similar approach has been used to calculate the prop- 
erties of several low-dimensional systems |2|, but so far 
no comparison with NMR-data of a direct realization of 
the AFHC has been done, especially in the vicinity of the 
QPT. 
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The NMR relaxation rate T x can be presented 
as T,- 1 '-2-2*2, 



{>yh%tf/2)Jdq[F*{q)S**{q,u N ) 
F x (q)S xx (q,ojN)], where j e and jn are the electronic 
and nuclear gyromagnctic ratios, respectively, u>n is the 
resonance frequency of nuclear spins, F z ' x (q) are the hy- 
perfine form-factors of nuclear spins, parallel and perpen- 
dicular to the external dc magnetic field B, S IJ ' v (q,ui]sr) 
{jjLv = z,x) are the components of the tensor of the dy- 
namical structure factor (DSF) of the AFHC, also par- 
allel and perpendicular to B. For the transverse com- 



ponents we have S xx = S vv , because of the rotational 
symmetry perpendicular to the field direction. Since 
ujn <C J/U,j e B, we use the limit uj^ — > 0. The ground 
state behavior of the DSF of the AFHC in the external 
magnetic field has been studied in [8j, where it was sup- 
posed that continua of low-energy gapless excitations of 
the AFHC (known as spinons) yield the main contribu- 
tion to the DSF. That conjecture, supported by numeri- 
cal calculations for small-size spin chains [8j], was later 
confirmed by the exact calculation, which used Bethe 
ansatz equations in the ground state of the AFHC [9(: 
About 75% of the DSF of the AFHC is determined by the 
two-spinon continuum. At B = the DSF is isotropic, 
and at w -> two points, q = and q = -k. con- 
tribute mostly to the relaxation rate. B ^ intro- 
duces anisotropy in the components of the DSF of the 
AFHC. For S / the main contribution to the com- 
ponents of the DSF of the spin chain is determined by 
the lower and upper boundaries of the two-spinon con- 
tinua [8| . The nonzero field causes the shift of the most 
important contribution for S zz {q,uj = 0) from q = n to 
q = 7r(l — 2m), where m is the magnetic moment per 
site of the AFHC, while the point q = remains impor- 
tant. For S xx (q,uj = 0) the point q — -k remains relevant, 
while instead of q = the main contribution comes from 

q = 27T771. 

The asymptotic behavior of the correlation functions 
of the AFHC can be calculated in the framework of the 
conformal field theory [10(. The low-energy states ap- 
proach zero (w = 0) at the vector q ~ Pp, where Pp is 
the Fermi momentum. For the longitudinal component 
of the DSF at q ~ Pf = we have 



nv 2 S* h * 



2Kak B T , 



(1) 



where K is the LL exponent, v is the Fermi velocity of a 
spinon, and a is the cut-off parameter of the theory. For 
q ~ Pp = 7r(l — 2m) we get 
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(2) 
y) is the beta function, 
and C\ is the field-dependent multiplier (ll|. At m = 
1/2 this contribution has to coincide with Eq. ([1]), which 
defines C\(B S ). The calculated longitudinal components 
of the DSF manifest weak dependencies on T and B, 
except in the vicinity of the QPT (at which v — > 0), 
where they show a strong growth linear in T. 

For the transverse component of the DSF at q ~ Pp = 
7T, we have 



vS xx 
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where 7' = 1/(AK). Near the saturation point B = B s 
the correlation amplitude goes to zero, while the value of 
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FIG. 1: The magnetic field (given in units of J/hje) and 
temperature (given in units of fes/J) behavior of the trans- 
verse staggered component of the DSF of the AFHC model, 
calculated using the conformal field theory. 



the correlation function at zero field in the ground state 
is approximately equal to 0.18 [11(. Hence, we can write 
the multiplier as C* 2 = 0.18/ J B 2 (l/4, 1/2) w 0.0065. This 
component yields the main contribution to the measured 
relaxation rate, see below, and the results of its calcu- 
lation are presented in Fig. 1. For q ~ Pp = 2irm we 
get 



.Wl + V^Wl,!. 



cos(27T7) 

-2Kak B T\ 2 ~ ) 



(4) 



where 7 = K — 1 + 1/(4 K). In the ground state nu- 
merical calculations |ll[ give the value of the correlation 
function at zero field ~ 0.03, which defines C%(B = 0). 
At B = the exponent for the AFHC is K = 1/2, [lfl| 
and, therefore, this component of the DSF, calculated in 
this approach, diverges at B = m = Pf = 0. However, 
that divergency is well-known to be nonphysical. It is 
easy to calculate the transverse homogeneous magnetic 
susceptibility for the Heisenberg spin system at q — 0; it 
is equal to m/Tij e B. For the AFHC the magnetic mo- 
ment is proportional to the field for small values of B, 
hence, in that region the transverse magnetic suscepti- 
bility coincides with the longitudinal one. The magnetic 
susceptibility is related to the DSF via the fluctuation- 
dissipation theorem. The correct magnetic field behavior 
of this transverse component of the DSF has to behave as 
the homogeneous longitudinal component for small val- 
ues of the field, and decay to zero at B — > B s , i.e. to coin- 
cide with the staggered transverse component there, be- 
cause at B = B s we have m = 1/2, and Pp = 2-nm = it. 
The velocity and the LL exponent of spinons are in- 
dependent. These dependencies can be obtained from the 
exact Bethe ansatz solution. Recently, a simple ansatz for 



the magnetic field behavior of the velocity v and exponent 
K, valid in the interval < B < B s , was proposed [12J : 
v = (7rJ/2 y [1 - (B/B 3 )] [l - (B/B s ) + (2h le B/7rJ)], 
K = //V4/ 2 - 3(n 7e S)2, / = ttJ[1 - (B/B s )} + 
h"f e B. The behavior of v and if, given by those 
expressions, agrees with the Bethe ansatz calcula- 
tions. Finally, marginal operators (from the renormal- 
ization group viewpoint) introduce logarithmic correc- 
tions to the asymptotic behavior of correlation func- 
tions of the AFHC in the conformal limit at low T 
[10| . Those corrections can be taken into account, 
see, e.g., [13(, which yields the additional multiplier 
^ln(24.27J/afc B T)/(27r) 3 / 2 to the right hand sides of 
Eqs. ©-©. 

During the last years a new approach for the calcu- 
lation of the critical exponents was developed 14j. It 
was pointed out that the low-energy dynamics of quan- 
tum chains is determined not only by Fermi points, but 
also by high-energy states of the system, i.e. the non- 
linearity of the dispersion relations was taken into ac- 
count. Contributions from those high-energy states were 
approximated as an interaction of the LL with an ef- 
fective impurity, which parameters are determined by 
quasi-momenta of excitations. The theory, which de- 
termines the renormalization of exponents of quantum 
chains in the presence of impurities, was presented in 
[lOj . Generalizing the approach of [14| . we conjecture 
that the B- and T-behavior of the DSF is determined by 
Eqs. ([I])-© with exponents, renormalized due to high- 
energy excitations. In the conformal field theory we re- 
place AM -> (AM - n imp ), AD ->• (AD - d imp ), where 
AM and AD are integers, determining the finite-size 
spectra of the chain, and rii mp and d imp are the param- 
eters of that effective impurity. The latter are defined 
as [H| n imp = ±(VK - 1), and d lmp = —(l/2y/K)m mp . 
The plus sign corresponds to the negative "valence" of 
the impurity, i.e. to the hole in the Fermi sea, and the 
minus sign corresponds to the positive "valence" , i.e. to 
the excitation above the Fermi sea. The behavior of 
the longitudinal homogeneous component of the DSF, 
Eq. ([I]), is, obviously, not renormalized. For the longitu- 
dinal staggcrcd-like component of the DSF, in Eq. $Z§ we 
need to replace the exponent K — > (1/2) — (l/2y/~K) + 
(1/4K) + (9K/16) - (3y/K/4) for a high-energy hole, and 
K ->• (1/2) - (l/2\/#) + (1/4A) + (if/4) - (\ r K/2) for 
a high-energy excitation. In Eq. ([3]) for the transverse 
staggered component of the DSF, we need to replace 
i -)■ (1/2) - (VK/2) + (if/4) for a high-energy hole, 
and i -> (1/2) + (1/if) - (1/2VA 7 ) + (A/4) - (VK/2) 
for a high-energy excitation. Finally, for the transverse 
"homogeneous" component, in Eq. (|4]) we need to re- 
place the exponent 7 -> (9A'/16) - (1/2) - (3y/K/4) 
for a high-energy hole, and 7 — > (1/K) + K — (1/2) + 
(\/if/2) — (l/2yK) for a high-energy excitation. How- 
ever, the T- and B-dcpendcncies of the components of 
the DSF, obtained within our conjecture (taking into ac- 
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FIG. 2: (Color online) Field dependence of Tr 1 at k B T/J = 
0.15 (error bars are within symbol size). The results of the 
effective field theory are given by the solid line for fields up to 
B a . Above B s an exponential function was fitted to the data 
(dashed blue line). 



count the high-energy excitation/hole of the AFHC), do 
not agree with the experimentally observed data. This 
can be explained as follows. For the calculation of the 
NMR relaxation rate one performs the integration with 
respect to quasi-momenta, where only a small interval 
near w = contributes mostly. The behavior of expo- 
nents in those intervals is not homogeneous: One expects 
the "traditional" LL exponent near the Fermi point and 
renormalized exponents close to the edges of the interval, 
with a smooth crossover between the exponents. The cal- 
culation of that crossover is a subtle point, not yet per- 
formed analytically. Thus, NMR experiments suggest, 
that in the interval of integration the main contribution 
comes from the region, where the "traditional" exponent 
is applicable. 

The NMR measurements were performed in two dif- 
ferent standard NMR setups, each with a supercon- 
ducting magnet, a 4 He temperature insert and a com- 
mercial/homebuild spectrometer. An inversion-recovery 
pulse sequence was used to measure the 13 C-nuclear re- 
laxation rate T^~ . In order to ensure a well defined 
comparability with theory, the hypcrfine form- factor F z , 
which relates T^ to S zz , was minimized. This was done 
by an orientation-dependent study of the NMR frequency 
shift and Tj -1 , determining the angle of l(B, b) = 50° in 
the b — c plane [5j . In this orientation the critical field is 
B s = {2J/hrf e ) = 14.6 T, which is a slight adjustment to 
the previous value, considering the anisotropic <?-factor 
from recently published ESR-results [15| . In this case we 
have Tf 1 = F s x {q)S* x + F' h x (q)S% x , which leaves only 
two form factors and the cut-off a as free parameters to 
fit our calculations to the experimental data. The results 
of this fit procedure are shown in Figs. [5] and G3 As for 
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FIG. 3: (Color online) Comparison of temperature-dependent 
experimental and theoretical Ti _1 -rates at different external 
fields. At low fields, an additional contribution from the ho- 
mogeneous DSF is included in the theoretical rates, see text. 



the independence shown in Fig. f2l the main contribu- 
tion to the Ti _1 rate comes from S xx . For fields larger 
than B s , a spin excitation gap opens linear with B — B s , 
leading to an exponential decay of the relaxation rate. 
We find an excellent agreement between our calculations 
and the NMR experiment for the whole region of fields, 
in particular near the QPT. 

The temperature dependence of T x ~ 1 at different fields 
is shown in Fig. [3J The values of F s x (q) and a obtained 
from the fit of the field dependence were used for the scal- 
ing of S xx at all B and T. At B close to B s , the experi- 
mental and theoretical NMR-rclaxation rates show a di- 
verging behavior as T — > 0. This singularity occurs when 
B, acting as the chemical potential for spinons, crosses 
the boundary of the dispersion relation (at this QPT v 
goes to zero), cf. Eq. ([3]). Again, the critical regime, i.e. 
T^\T) at 12.8T and 13.8T, can be fully described by 
F s x (q)S xx . Fitting the experimental rates with a power- 
law decay for T < 2J/Ti-y e , in the region 0.7 J < T < 2 J 
we find the exponents to be slightly smaller than the 
theoretically calculated values , cf. Fig. [3] Notice, that 
one expects the accuracy of field-theoretical calculations 
to significantly decrease for T > J/ks- The agreement 
with experiment is, nevertheless, surprisingly good for 
the whole region of temperatures measured, with T x ~ be- 



coming almost T-independent at high temperatures. At 
low fields, the experimental rates show an approximately 
linear T-dependence up to T ~ 2J/fcs- This behavior 
is not reflected by F s x (q)S xx . As described above, for 
low fields S% x has to be equivalent to S^ z . The fit of 
1/Tj = F' s x (q)S xx + F' h z (q)Sf l z for T< 2J/Ti le , with F h z 
as a free parameter gives much better agreement with 
the experimental data in this low-T region. The ratio of 
F h x (q) at 2T and 6T is ~ 5/2, explainable by the shifting 
of the Fermi point in q as a function of B and the asym- 
metric position of the 13 C probe between two magnetic 
Cu 2+ ions. 

To summarize, based on the conformal field theory, we 
presented a comprehensive calculation of both the trans- 
verse and the longitudinal DSF of the AFHC model in 
the LL regime. Considering the hypcrfine coupling as a 
free parameter, the comparison to NMR results on the 
AFHC system CuPzN was performed for temperatures 
0.19 < {k B T)/J < 4.5 and fields up to B ~ 2.2J/fv^ e . 
The spin dynamics close to the quantum critical point 
can be extraordinary well described by the staggered part 
of the transverse DSF. At low fields contributions from 
the homogeneous part of the DSF have to be taken into 
account. We thank T. Giamarchi for suggesting us to 
perform this study, L.I. Glazman for the discussions on 
the range of applicability of exponents, and W. Brenig 
for the discussion on the QMC. This work was supported 
by the DFG through Grant No. KL1086/8-1 of FOR 912. 
AAZ acknowledges the financial support by the DFG via 
the Mcrcator program. 
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